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SUMMARY 


A two-degree-of-freedom system, in which rolling oscillations are 
forced, is described. Details of the system, the theory, and the method 
of operation are discussed. It is shown that, although the system is 
characterized by nonlinear equations of motion, Linearization of the 
equations by assuming small perturbations and constant coefficients yields 
sufficiently accurate results. 


The accuracy of the system was first investigated in tests under 
idealized conditions in which the aerodynamic derivatives were simulated 
by the action of gyroscopes and magnetic dampers. Later, the system was 
investigated in a high-speed wind tunnel using a simple model for which the 
aerodynamic derivatives could be estimated. 


Lhe results of the tests showed that the quantities relating to the 
primary mode of operation of the system, that is, the rolling velocity 
derivatives, could be obtained satisfactorily. The same is true of the 
directional stability and the damping-in-yaw derivatives. Results obtained 
from the data-reduction equations for the rolling moment due to yawing 
velocity and due to sideslip angle were unreliable, however, and prevented 
the evaluation of these derivatives. 


* 


INTRODUCTION 


The increasing emphasis placed on studies of the dynamic stability 
of aircraft has intensified the need for satisfactory methods for predict-~- 
ing the dynamic stability derivatives. Most theoretical methods evolved 
to date are subject to limitations imposed to reduce the mathematical 
task to reasonable proportions. In order to obtain experimental checks 
of these methods and to provide sources of information in those areas in 
whick theoretical methods are not yet available, there exists a need for 
experimental research on the stability derivatives. 


The experimental study of dynamic stability derivatives in wind tun- 
nels is, of course, not new. Many systems have been evolved for this 


purpose. References 1, 2, and 3, for instance, discuss early British 
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systems in which the relative motion of the model and the air stream is 
generated by model motion. A dynamic system recently developed at Ames 
Aeronautical Laboratory is described in reference 4; again the relative 
motion of the model and the air stream 1s generated by model motion. 

The rolling- and curved-flow wind tunnels now in use at the Langley Aero- 
nautical Laboratory of the NACA, for instance (refs. 5 and 6), are examples 
of other systems in which the model is held rigidly and the air flow gen- 
erates the relative steady yawing and rolling motion.. 


In the present case, a system in which the model executes forced 
oscillations with two degrees of freedom was developed. The details of 
its design and development are published herewith since they may prove 
useful to others who desire to develop dynamic systems. 


SYMBOLS AND COEFFICIENTS 


The aerodynamic coefficients defined herein are referred to the 
system of body axes originating at the model center of gravity. The 
symbols and coefficients are defined as follows: 


b body length, ft 
a 1 2V 
2V 
“Ly Y \ asb? 
l 
“2g P (+ 
: 2V 
“28 Lé (Æ) 
. 2V 
“np "p ~*~ 
” 2V 
Cny Ny > 
a 
Cng Ng > 
. fev 
ng Ng at) 


Í frequency, ceps 


Tl 
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Hy Hy, Hy, 
hy 3 hy Nz 


Ix 


moment of momentum about space-fixed axes 

moment of momentum about body axes 

moment of inertia about x axis, ft-lb sec? 

moment of inertia about z axis, ft-ib sec? 

product of inertia with respect to x and z axes, ft-lb sec? 
product of inertia with respect to y and z axes, ft-lb sec? 
rolling moment, ft-lb 


absolute value of spring constant of roll flexure pivots, 
ft-lb/radian 


rate of change of aerodynamic rolling moment with rolling 
angular velocity, ft-lb/radian/seec 


equivalent viscous damping factor due to internal friction in 
roll flexure pivots, ft- lb /radian/sec 


rate of change of aerodynamic roLling moment with yawing 
angular velocity, ft-lb/radian/sec 


rate of change of rolling moment with sideslip angle, 
ft-lb /radian 


rate of change of rolling moment per unit rate of change of 
Sideslip angle, ft- 1b/radian/sec 


pitching moment, ft-lb 
applied moments about space-fixed axes, ft-lb 
yawing moment, ft-lb 


rate of change of aerodynamic yawing moment with rolling 
angular velocity, ft-lb/radian/sec 


absolute value of spring constant of yaw flexure pivots, 
ft-1b/radian 


rate of change of aerodynamic yawing moment with yawing 
angular velocity, ft-lb/radian/sec 


equivalent viscous damping factor due to internal friction in 
yaw flexure pivots, ft-lb/radian/sec 


rate of change of yarns moment with sideslip angle, 
ft- lb/rađian 


Neg . atts 
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rate of change of yawing moment per unit rate of change of 
sideslip angle, ft-lb radian/sec 
rolling angular velocity, radians/sec 


dynamic pressure, ib/sq ft, sev? 


maximum body radius, ft 

yawing angular velocity, radians/sec 

maximum cross-sectional area of body, sq ft 
roll input torque, ft-lb 

time, sec 

free-stream velocity, ft/sec = 

Space-fixed axes 

axes fixed with respect to the model (body axes) 
angle of attack, radians 


angle between free-stream velocity and longitudinal axis of 
support housing i | 


angle of sideslip, radians 7 
angle of roll, radians 
angle of yaw, radians 


phase angle between roll and yaw angles, positive for yaw 
leading roll, deg : 


phase angle between roll angle and roll-input torque, positive 
for input torque leading roll, deg . 


angular frequency, 2xf, radians/sec 
mass density of air, slugs/cu TE 
amplitude of oscillation or out-of-phase component 


initial condition or in-phase component 


al) 


dt 


Fi 
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ANALYSIS AND DESCRIPTION 


Coordinate System 


The coordinate system for which the equations of motion will be 
written is shown in figure l. The x,y,z axes are fixed with respect 
to the model and are therefore body axes. The X,Y,Z axes are alined 
with the equilibrium position of the model and remain fixed in space. 


In order to illustrate the possible modes of motion of the model, 

a schematic drawing of the mechanical apparatus is shown in figure 2(a). 
As shown in this diagram, the main shaft (torque tube) of the apparatus 
is driven in a constant-amplitude oscillatory rotation about the X axis 
Which is transmitted directly to the model. The axis of rotation of the 
torque tube is defined by the roll flexure pivots and ball bearing which 
support the tube within the support housing. Attachment of the model to 
the torque tube is through the yaw flexure pivots which allow the model 
a second mode of motion, that of angular rotation about the z axis. 


Equations of Motion 


In accordance with the Newtonian laws of motion, the time rate of 
change of the moment of momentum about axes fixed in space is equal to 
the summation of the externally applied moments about those axes: 


> ma 

i as dt 
dity (1) 

5 pee 3 

E at 

cHz, 

SMe = — 

Z at 


Since the model is rigidly constrained in pitch, the second of equa- 
tions (1) need not be considered. 


The time rate of change of the moments of momentum about fixed- 
space axes in terms of the rate of change about the body axis system may 
be written (see, i.e., ref. T): 


where 


hy = 


The applied moments and inertial 
of the body axis variables 


OL 


XN 


il 
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a ee 
(2) 
Z e 
——— 
at y F 


Ply + Vz 
-PI xy Pa YIyz (3) 
plz $ YIz 


reactions can then be expressed in terms 


Tx + Ixg¥ + Inyo + Iy% 


we a a> ae 
Ig¥ + Ixz? - Ixy - Tyz¥ 


Substituting the externally applied moments for the left-hand side of 


equations (4) gives 


T cos Y + Los? cos V + Lof cos ¥ 


- Lop + LYY + LAB + Leb 


20-4 igt d (Txyoh + ta (5) 


Wye + Nyy + NOP - Npk + NAB + NgB = Te + IxzỌ + (-Iyy - yh) 


A geometrical relationship exists between the angle of sildeslip B 
and the roll and yaw angles @ and y as follows: Let the angle between 


the free-stream velocity and the 


š$ axis be Qile 
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From sketch (a), 
w = V sina Vi = V cos ay 


and from sketch (b), 


v =W sin Ọ= V sin aı sin Ọ 


From sketch (c), 
Sketch (a) 


vı = V cos y = V sin a, sin Ọ cos y 


Vo = Vi sin ¥ = V cos a, sin y 


Tke angle of sideslip is then given by 


TD 

UI 

ve 

Je 

rs 

? 

bh 
o> 
} 
+ 
< 

N 
2 a i 


(6) > 


B = sin™?(sin a, sin Ọ cos ¥ - cos aq sin y) Sketch (b) 


The true angle of attack is also a 
function of the roll angle: 


sinl (sin az cos ¢) (T) 


Substitution of equation (6) and its 
time rate of change in equations (5) 
yields the equations of motion, and 
it can be seen that nonlinear terms 

Sketch 
Will appear both in the aerodynamic etch (c) 
moments and the inertial reactions. 


in order that they be of practical value for use with experimental 
wind-tunnel data, linearization of the equations is necessary. If the 
approximation is made that 
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sin 0 = © sin y= y 
cos Ọ = cos y = l Po 


then the angle of attack becomes a = a1, and the sideslip angle is given 
by 


B = @sina - y cosa (8) 


In general, the angle of attack is sufficiently small so that equation (8) 
may be written as 


p=ap-y. (9) 


and the time rate of change of sideslip angle is given by 


Bb =ap - ¥ (10) 


For large angles of attack for which these approximations are not valid, 
equation (8) and its time derivative may be used. r 


If, in addition to the approximations made in equations (9) and (10), 


the second-order inertial terms in equations (5) are neglected, the equa- ia 
tions of motion may be written 


lp tia 


T i (Lo, + Lb + abs) ® - (Lo - aLg)p + (Ly - L3)¥ - Ley 
(11) 
Tew + Tue 


(Nọ + aNg)ó + aNgp + (Ny, + My - Nay - (Ny + gly 


The question immediately arises as to the accuracy with which equa- 
tions (11) represent the physical system. In order .to answer this question, 
transient solutions of both the nonlinear and linearized equations of 
motion for a pulse-type roll disturbance T _were obtained on the Reeves 
Electronic Analog Computer for conditions expected to be encountered in 
wind-tunnel operation. Comparison of the solutions showed that the only 
discernible effect of linearization was a negligible difference in the © 
roll and yaw accelerations. From this result it was concluded that the 
use of equations (11) would be a justifiable approximation, and that their 
use would yield data of sufficient accuracy. Further experimental veri- 
fication of the validity of these equations is given in a later section. 
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Solution of Equations of Motion 


In the following development the solution of the equations of motion 
will be given and arranged in a manner which will permit the determination 
of the desired aerodynamic parameters in terms of measured quantities. 


Yaw equation.~ The second of equations (11) may be rearranged and 
written 


IZY - (Np, + Nọ - NAW + (Ny + Np) ¥ = - Ix + (Nọ + aNG)$ + aNgM (12) 


Equation (12) is a linear second-order differential equation, the solution 
for which is well known. The roll angle ọ is forced to vary harmonically 
with time; Ọ= 9% sin wt. Then the solution, in general, consists of twa 
parts: first, the complementary function, uniquely determined by the 
initial yaw angle ýy and the initial yawing velocity w; and second, the 
particular integral, sometimes called the steady-state or frequency 
response, which is independent of the initial conditions, and is deter- 
mined only by the applied disturbance, the right-hand side of equation (12). 


The amplitude ratio and phase angle given by the particular integral 
are, 


Yo (aNg + Izt) + (Nọ + aNg) oF 133) 
o (Ny + Ng - Igw?)? + (Ny, + My - NB) ?w? 
No + aNg )w Ve + Ny - Nå 
E = tan7+ in li 4 tan +t ys + Ny - N)0 (14) 
aNg + Ix2W (My + Ng) - Iæ? 


Equation (14) may be manipulated to give 


w (Ng + ong) | (My a Nee = T,w)cos € + (My. + Ny - Ng)w sin e] 


aNg + AA = 3 
(My + Ng - iol) Sin € -= (My, + My - Ng)w cos € 
(Nb + aNg) (aNg + Ixzw~) R? + Ng - Izw=)sin ¢€ - (Nis + Nir - NÉ Jw cos e] 
! gjo = — MMM I 


(Ny + Np - T,W=\)eos e + (My. + My - Np )w sin € 
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Substitution of these equations into equation (13) leads to two funda- 
mental relationships: 


Yo (Nọ + aNg )w , 
"oO _ 7 15 
Po (Ny + Ng - I,w-)sin € - (Ny. + Ny - NÅ Jw cos € 

Vo aNg + Tyg 

a. = (16) 


Po T (iy + Np - TpuA)eos e + (ij, + Ny ~ Ng)w sin € 


The geometrical meaning of equations (15) and (16) is apparent from the 
concept of the rotating vector diagram (see, e.g., ref. 8). Since 9 

is a harmonically varying function 9 = % sin wt, the right-hand side — 
of equation (12) represents a harmonically varying yawing moment No 

with components proportional to the roll position, velocity, and accelera- 
tion. The vector diagram is shown in sketch (d). The magnitude of the 


resukting yawing moment is 
(Agt ang) Pr w 


No 
£ 2\2 ° : 2 
baT ake a = 


_, No + aNg Jw 


utr 
iI 
cr 


aNg + Ingzw 


The particular integral shows that y is also a harmonic function, 
and equation (12) states that the sum of the moments R, resulting from 
(N; +My - Ng] Yoo the motion Y must be equal and 

i opposite to No. The vector diagram 
Lha? is shown in sketch (e). The magnitude 
(Ny +Ng) Yo Soa is of the resulting yawing moment No | 

(which is equal to R) is 


and the angle p is given by 


(Ny, + Ny - Ng) 


Ee E E 
5 Ny + Ng - Izo? 
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The. amplituđe ratio and phase angle given by equations (13) and (14) may 
then be written 


Yoa No/Vo 
o = No /¥o 
= re Se 


or 
(Ny +N g -Ng) faw j PATI. 

Equation (12) may now be represented 
by combining the vector diagrams shown 
in sketches (d) and (e}. It can be 
seen that equations (15) and (16) 
represent the summations of the out- 
of-phase and in-phase components of (My +No) te 
the vector diagram of sketch (f). Sketch (f) 


The quantities in equations (15) and (16) measured during operation 
of the dynamic apparatus are o; Yo, w, and c€. At any discrete frequency 
these are insufficient data to solve the equations; however, use of the 
apparatus as a free-oscillation system provides additional data from the 
complementary solution of equation (12): 


(Ny. + Ny - Ng)t (My. + My - Np) 
¥ = \ exp [ERLE H feos wat = oe sin vat | a 


ut sin vat} (17) 


where the angular frequency of oscillation is given by 


pty + Ne E e ae ary) a 
ve Iz Plz 


For wind-off conditions, equations (17) and (18) may be written 


Nut Ny Y; 
y = lex mad e) | va (cos wnt - a sin ont) + TA sin ont (19) 
Z Zn 


where 
_ | _ (ey | (20) 
PE A te 2T 
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The contribution of the system damping to the oscillation frequency is 
negligible, and equation (20) may be written 


Ny 
oT 


The moment of inertia may then be evaluated: 


Ny | 


The system damping may be determined as 


š aa 
Ny = (22) 
5 Wind off 


where Ve and Y are ordinates of the envelope curve a time interval T 
apart; simi ariy, the aerodynamic damping becomes 


el, i 
Ny - NA =— in s Ne (23) 
ý : i Vt/wind on Ys 


and the directional stability parameter Ng may be determined from 
equation (18): 


(My ia Ny = na)” 
¢.om 


Ng = 4n*I,fa* i, 


- Ny (24) 


Sufficient data are now available to solve equations (15) and (16) for 
wind~on conditions: Sa 


yn (Ny + Ng - I,w-)sin e - (Ny, + Nọ - Nj)w cos e€ 
Ng + aNg = a gee ee = (25) 


Vo (Ny + Ng - Iew-)eos e + (Nys + a - Ng )w sine aNg 


xz = 
= w= We 


ill 
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It should be noted that the product of inertia fy, may be obtained 


for wind-off conditions by elimination of the damping term in equa- 
tions (15) and (16) to give 


o Yo (My = to) 
2 Po w=eos € 


However, it was found that, for winged models at angle of attack, the 
Slight elastic deformations of the model and support structure were suf- 
ficient to necessitate the use of equation (26). 


Roll equation.- The first of equations (11) may be rearranged and 
written as follows: 


I% - (Los + Lo + ag) g + (Le. - abe) P + Tay É (Ly - Lg) +igf =T (27) 


Before proceeding with the total equation (27), we will consider the case 
of single-degree-of-freedom roll motion, defined by equation (27) when y 
and its derivatives are identically zero. We proceed directly from the 
vector diagram given in sketch (g). 

Summation of the out-of-phase and in- (Lg tlgtalg) foe 


phase vector components gives ü 


A 
e P E LEE s 
(Lo, + Ig + alg)Pow + To sin A = O (28) Lap- w! algfe Iseo? i 


5 
Izbu? - (Lo - alg) Po + To cos A=0 (29) ee: 


For wind-off conditions, equations (28) and (29) may be solved for 


. To sin A 
Ss eee (30) 
-To cos A + La 
Ty = = (31) 
Dow 
For wind-on conditions, 

lo Sin A | 

ë > S a ee eh odis at g a a 

Lo + ala Í Pou Los (3 


1 To cos A + IyPo 
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Equation (33) becomes indeterminate at zero angle of attack; however, 

another expression for Lg may be 

obtained from the more general equa- 3 
tion (27) which we now consider. 

The vector diagram representing 
equation (27) is shown in sketch (h). 
Summation of the out-of-phase and 2 
in-phase vector components gives 


We 
A Sketch (h) 


To sin A + TxzVow“sin € + (Lo. + Leo + ALG) Pow + 


(Lý - Lg)¥o@ cos € - Lg¥o sine =0 (34) 


To cos A + (aLg - Lo) 9g + IxPow? + Iyz¥ow cos € - : 


(Ly - Lg)¥,ow sin e - Lg¥, cos € = 0 (35) 


Multiplying equation (34) by sin e, equation (35) by cos e, and adding 
and solving for Leg gives 


l 5 
ap. cos € - Yo ~ g - - we m 
Lg a0, cos E -= Yo o Lyx & )®o cos € To cos(A €) IxzYo 


(Las + Lọ + abs) Pow sin e| (36) 
where the product of inertia I,, is given by equation (26). Multiplying 


equation (34) by cos e, equation (35) by sin ec, and subtracting and 
solving for (Ly - LA) gives i. : 


Ly - Lg x (cg - Lo + Ixu*)@, sin e - (LO, + Lo + aLg) pow cos € -~ 
Ty sin(A - e) | (37) 


Description of Apparatus i 


A physical description of the drive apparatus and its various com- 
ponents is presented below. The next section deals with the functions of 
the instrumentation. | S ; 
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In discussing the drive apparatus, it is convenient to refer to the 
pictorial representation of the system presented in figure 2{b). The 
model is allowed two angular degrees of freedom (rolling and yawing), 
subjected to constant-amplitude forced oscillations in roll, and is free 
to respond in yaw. It is rigidly restrained in pitch (at a constant 
angle of attack from 0° to 15°), and all externally applied pitching 
moments, lift, drag, and side forces are transmitted directly to the 
wind-tunnel structure. The apparatus is powered by a three-phase induc- 
tion motor mounted vertically over the support-strut housing at the rear 
of the apparatus. Rotary motion of the motor shaft is converted to an 
approximate harmonic motion of the main drive shaft (torque tube) by means 
of the eccentric block and drive flexure indicated in the figure. An 
enlarged drawing of these components is shown in sketch (i) to aid in 
visualizing the principle involved. 

Amplitude of the oscillatory rolling | 

motion is determined by the amount of 
eccentricity in the eccentric block, and 
with the present system is 40. ‘The 
torque tube, which transmits the oscil- 
latory motion to the model, is supported 
within the stationary support sting by 
the roll flexure pivots at the forward 
end and by a ball bearing at the rear. 
Rigidly attached to the front end of the 
torque tube are the yaw flexure pivots 

on which the model is mounted. The only 
contact between the model and the appa- 
ratus is through the mounting at the for- 
ward end of the yaw flexures; sufficient 
clearance exists between model body and support sting to allow the model 
to yaw approximately 2° about the flexure axis. The maximum operating 
frequency of the system is approximately 20 cycles per second. 


Sketch (i) 


To obtain accurate results with the system, it is necessary that 
rapid variations in oscillation frequency of the apparatus be eliminated. 
This purpose is served by the small flywheel shown mounted over the drive 
motor in figure 2(b). 


The use of flexure pivots in the apparatus resulted primarily from 
the need for avoiding the often undesirable damping characteristics of 
ball bearings. The flexures used were designed in accordance with infor- 
mation presented in reference 9, and were found to be entirely satisfactory 
from the standpoint of defining a rigid axis of rotation. 


Characteristics of flexure pivots.- The flexure pivots are designated 


as compensating type in reference 9 because, theoretically, the spring 
constant is essentially independent of load. The desirability of this 
feature in the present application is obvious, and the experimental veri- 
fication is shown in figure 3 for the case of the yaw flexure pivots sub- 
jected to lift loads of 0, 100, and 200 pounds. The invariance of the 
roli-flexure-pivot spring constant was equally good. 
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The internal damping in the yaw flexure pivots was obtained by means 
of the well-known oscillation-decay method. It was found that the rate 
of decay was not logarithmic but, instead, wes a function of both ampli- 
tude and frequency. However, equation (22) was applied to the data in 
order to obtain an equivalent viscous damping factor Mg It was found 
that, as shown in figure 4, the data could be represented by an equation 
of the form 


š a + bY 


Ny. = F E ; (38) 


Justification for this type of expression can be seen by converting to 
the form given in reference 10 where the internal damping is presented 
in terms of a "damping capacity" Aw, the damping work per cycle per unit 
volume. This quantity was found to be independent of frequency and to 
vary as the third power of the stress amplitude. Consider the case of 
constant-amplitude oscillation, ¥ = yo sin wt. Then the damping capacity 
may be written (see, e.g., ref. T) with the flexure-pivot cross-sectional 
area A and length h, | 


2x ; 
x 2,2 — oN 2 
ee ete a ie cos“wt dt = ee 
Ah iS Ah 


or; 


(a + bio) 207h* 


a (39) 


Although Aw was found to vary in reference 10 primarily as the third 
power of the stress amplitude (which is directly proportional to the 
oscillation amplitude), examination of the results therein indicates that 
at the lower amplitudes the data may be approaching a second power vari- 
ation as indicated by equation (39). Likewise, the damping capacity 
given by equation (39) is independent of frequency, and it is felt, there- 
fore, that the use of the equivalent viscous damping factor represented 
by equation (38) is justified. The data presented in figure 4 also show 
internal damping for the case of the yaw 
flexure pivots subjected to lift loads of 
O, 100, and 200 pounds. The application 
of an external_load which produces no addi- 
tional damping is essentially impossible. 
However, it was found that the additional 
damping of a lift load simulated as shown 
in sketch (j) was negligible, and the data 
were obtained in this manner. Again it can 
weights be seen that there was essentially no effect 
| of Lift Load. 
Sketch (5) 


N 
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No successful method was devised for the application of a lift load 
which provided no additional damping about the roll axis. Consequently, 
the internal damping of the roll flexure pivots was assumed to be invariant 
with lift load, an assumption justified by the yaw-flexure results, and 
the fact that any small change that might occur would be a negligible 
percent of the aerodynamic roll damping. The equivalent viscous damping 
factor Lg, for the roll flexure pivots was obtained by operating the 
apparatus under wind-off, tunnel-evacuated conditions, and it was found 
that these results could also be represented by an equation of the form 
of equation (38). 


The values of Id, and Nj, were found to be approximately 20 percent 
and 12 percents respectively, of the aerodynamic values obtained in the 
wind-tunnel investigation to be described later in the report. It should 
be noted, however, that the aerodynamic roll damping measured for the 
body-fin model used in the wind-tunnel test is much smaller than would be 
the case for a complete airplane model. 


Tostrumentation 


Three quantities are measured by the system: angular roll and yaw 
positions of the model, and the roll-input torque applied to the model 
to maintain the forced oscillations in roll. These quantities are measured 
by means of resistance-type strain gages mounted on suitable beams, the 
locations of which are indicated in figure 2(b). The amplitudes and 
phase angles of the oscillating strain-gage signals are obtained by pass- 
ing each of the signals through a sine-cosine resolver located at the 
drive motor and gear driven by the motor shaft (see fig. 2{b)). The two 
output signals of each resolver are then directly proportional to the 
in-phase and out-of-phase components of the respective strain-gage signals 
referenced to the resolver rotor position. By the vectorial summations 
of these quantities the required amplitudes of oscillation are obtained, 
as well as the phase relationships necessary to establish the roll-input- 
torque and yaw-position phase angles with respect to the model roll posi- 
tion. Oscillation frequency is measured by means of an electronic counter 
driven by a portion of the oscillating roll-position strain~-gage Signal. 


Circuits.- It was noted in the preceding section that it is necessary 
to know the amplitudes of the roll angle, roll-input torque, and yaw angle, 
and the phase angles of the roll-input torque and yaw angle with respect 
to the roll angle. The desired quantities were measured by means of three 
identical electronic circuits, the essential features of which are shown 
in sketch (k). 


18 


d.c. voltage 


reference 


strain 400 ~w 


reference 


[ \ tal il | | ic 
g Aa | | hi 
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stotic deflection 


e^ 


galvan- 
ometer 
galyon- 
ometer 


stotic defiection 


Sketch (k) 


The principle of operation of the circuits is as follows: 


A direct 


voltage is impressed on the strain-gage bridge which is to provide a 


measurement of the desired quantity. 


The resulting output is an alternat- 


ing voltage whose frequency is that at which the system is operating and 
whose magnitude is directly proportional to the amplitude of the desired 
quantity. This output is used to modulate a 400-cycle wave; the result 
is amplified and then passed through the resolver which reduces it to 
components in phase and out of phase with the position of the resolver 
rotor. These components are then demodulated and used to drive the gal- 


vanometers as shown in the sketch. 


Resolvers.- The resolvers utilized in the apparatus were of the 
standard induction type, commercially available. Two coils are accurately 
positioned at 90° with respect to each other on the stator and, with a 
l-to-l transformer ratio, produce two harmonic signals also phased 909 
when the rotor is turned at a constant angular speed. 
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The action of the resolver 
may be explaimed as follows: 
Assume the quantity to be measured, 
as represented by the oscillating 
strain-gage signal (hereafter 
referred to as the basic signal), 
has an amplitude of oscillation, A, 
and phase angle, ©, with respect 
to the resolver rotor. This sig- 
nal is then defined as A sin(wt+5). 
The resolver may be considered to 
provide two reference signals of 
unit amplitude, sin wt, cos Wt as 
shown in sketch (2). Transmission 
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of the basic signal through the resolver effectively results in multipli- 
cation of the basic signal with the two reference signals: 


Vi = A sin(wt + 5)sin wt 


= = cos 5 - Š cos(2wt + 8) (40) 
Vo =A sin(wt + &)cos wt 
5 = sin 5 + = sin(2wt + 5) (41) 


Thus the amplitudes of the resolver output signals are one half that 
of the basic signal, the frequencies are twice the frequency of operation 
of the system, and each of the resultant signals has been displaced, the 
first by an amount equal to one half the product of the basic signal and 
the cosine of the phase angle, the second by one half the product of the 
basic signal and the sine of the phase angle. 


Equations (40) and (41), which represent the components of the basic 
strain-gage signal in phase and out of phase with the resolver rotor 
position, also represent the time-varying voltages used to drive the 
galvanometers as shown in sketch (k). 


Galvanometers.~ A galvanometer element is essentially a single- 
degree-of-freedom system, and its amplitude and phase-angle responses to 
a harmonically varying forcing function are well known, If the frequency 
given in equations (40) and (41) is much greater than the undemped natural 
frequency of the galvanometer element, the response of the element to that 
portion of the signal will be essentially zero. The first term on the 
right-hand side of each of equations (40) and (41) represents a voltage 
which is independent of time, however, and the galvanometer is capable 
of responding with a deflection G which is directly proportional to 
the static signal. The results, as shown in sketch (k), are static indi- 
cations of the in-phase and out-of-phase components of the harmonically 
varying quantity. 


Application of the proper calibration constants Co and Ci now per- 
mits the computation of the amplitude and phase angle of the desired 
quantity, the amplitude being equal to the square root of the sum of the 
squares of the in-phase and out-of-phase components, and the tangent of 
the phase angle being given by the ratio of the out-of-phase to in-phase 


components: 
A = 2 /03"G4* + Coto? 


(42) 
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The phase angle thus obtained are those between the resolver rotors 
and each of the measured quantities; roll positior, yaw position, and 
roll-input torque. It will be recalled, however, that the desired phase 
angles are those between the yaw position and rell positions, and between 
the roll-input torque and roll position. It is clear that these quantities 
ere immediately obtainable from the above information if the relative phase 
angles of the three resolvers with respect to each other are known. Since 
it is convenient to the operation of the system to maintain these relative 
phase differences at small values, the resolvers are first alined with each 
other as accurately as possible. A calibration is then made from which the 
remaining phase angles are obtained. ey 


Effects of harmonics.- When the principle of operation of the resolver 
circuits was illustrated, it was assumed that the basic signal was a pure 
harmonic function. In general, the quantities to be measured are not pure 
harmonic functions, but may be distorted due to extraneous disturbances 
from support vibration and other sources, one of which wes the slightly 
impure harmonic motion of the model produced by the eccentric block and 
drive flexure in converting the rotary motion of the motor into an oscil- 
latory motion of the torque tube. In addition, the voltage output of the 
strain-gege bridge may have a mean value other than zero due to initial 
unbalance of the bridge. The voltage output may then be represented by 
the Fourier series 


Ay + . (a, cos nwt + bn sin nwt) 


=1 
or, in.a@ more convenient form, 


co 


hat J An sin(nwt + ôn) (43) 
neo N21 sae l l 
An = Ee: + bn? 
1 D 
Syn = tan s 5 


In a similar manner, the resolver reference signals may not be pure har- 
monic functions; if, in addition, the two reference signals are not phased 
exactly 90° with respect to each other, but have some other phase angle 7; 
the Fourier series representing the two signals may be written 


is.) 


J Bm sin mutt (4h) 


m=1 
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(s) 


a By Sin(mwt + 7) (45) 


m=} 


Then the in-phase and out-of-phase voltages used to drive the galvanometers 
will be given by the products of equations (43) and (44), and (43) and 
(45), respectively: 


a co 0 | 
Vi = Ay - Bm Sin mwt + » Y Anba sin(nwt + 5y)sin mut 
m=1 


n=1 m=1 
ve) ; o co 
= Ag ) Bq sin mut + > X = feos| (2 - m)wt + On = 
m=1 n=1 M=1 
cos |(n + m)wt + òa} (46) 


Ag X Bm sin(mwt + y) + J » AnBm Sin(nwt + Sy)sin(mut + 7) 


Vo = 
m= n=1 m=1 
= AnBn 
= ho) By sin(mut + 7) + ee cos|(n - mut + 8a = 7) Š 
m=i n=1 m= 
cos| (n + m)wt + Sy + ri} (47) 


Examination of equations (46) and (47) reveals that time-independent 
voltages occur only When n = m, that is, for products of equal harmonics. 
Then the effective voltages, those to which the galvanometers will respond, 
may be written 


5 
Vi = » = AnBn (48) 
n=1 
= cos{(5y - 
vo = ) SRA? amn (49) 
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These equations show that the in-phase and out-of-phase components of 
each measured quantity may be in error by an amount equal to the sum of 
the products of equal harmonics contained in the quantity to be measured -o 
and in the resolver signals. If either of the signals has no harmonic l 
content, but consists only of the fundamental (n = 1}, no error is caused = 
irrespective of the distortion present in the other. Under normal wind- 

tunnel operating conditions it would be expected that, due to previously 

mentioned disturbances, some distortion would be present in the basic - 
signalis, and it is necessary that the resolver signals be free of all _ 
comparable harmonics in order to eliminate this possible source of error. 
For this reason the resolvers were driven by 90-tooth gears. A remaining 
possible source of low-frequency harmonic content in the resolver signals 
was slip ring "noise"; a visual examination of the resolver reference 
signals by means of an oscilloscope revealed no distortion, however. 

One other possible error due to the resoivers is indicated in equa- 

tion (49): the phasing of the resolver reference signals with respect 

to each other by some angle y other than 90°. For the resolvers used 
in the present investigation the angle y was accurate to 90°, +1 /2°, 

and the error from this source was, therefore, negligible. Now, since the 
fundamental amplitude, Bı; of the resolver reference signals is unity, 
equations (48) and (49) simplify to. 


A 

w= FA cos ô; a 
A E 

Vo = 5 sin ô, 7 


The galvanometer deflections G; and G, are directly related to the in- 
phase and out-of-phase voltages by the appropriate calibration constants 
Ci and Co. The amplitude and phase angle are given by 


Ay =2,/ Cate + C Ga 


Equations (50) are identical with equations (42) which were developed 
for the case of no distortion in the quantities to be measured. It can be 
seen that the system, in effect, filters out the unwanted harmonics end 
provides a measure: only of the desired response. 
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Frequency characteristics.- Mention was made of the alinement and 
calibration necessary to determine the final phase relationships of the 
resolvers. The method of alining the resolvers is as follows: While 
the oscillation frequency is maintained at some constant value, each of 
the three circuits is, in turn, connected to the output of the roll- 
position gage. In each case the appropriate resolver housing is then 
rotated until the out-of-phase component indicated on the galvanometer 
is approximately zero and is locked in this position. In. conjunction with 
these operations a check was also obtained on the frequency characteristics 
of the measuring circuits. It should be recalled that the phase angles 
necessary for data reduction are measured with respect to the roll posi- 
tion. Then it is clear that the absolute value and frequency dependence 
of the phase angles indicated by the resolver circuits is immaterial if 
these characteristics are identical for each circuit. The procedure was 
as follows: After alining the resolvers, each of the circuits was, in 
turn, connected to the output of the roll-position gage, and the ampli- 
tude and phase angle were computed for discrete frequencies covering the 
available range of operation. It was found that the phase shift varied 
Slightly for the three circuits, the phase angle differences varying in 
a manner linear with frequency. The phase shift difference was then used 
as a calibration for correcting the phase angles obtained during testing. 


Operation of the System 


The operational procedures followed in conducting wind-tunnel inves- 
tigations with the testing apparatus are described briefly below. It 
should be noted that the procedures are essentially identical when the 
apparatus is used for calibration tests such as are described in a later 
section. 


(1) With the model rigidly restrained in yaw, measurements of roll 
position 9, roll-input torque T, and the roll-input~torque phase angle 
A are obtained for single-degree+of-freedom roll oscillations, at dis- 
crete frequencies throughout the available range, by means of the resolver 
circuits previously described. Data thus obtained, for both wind-off 
tunnel-evacuated, and wind-on conditions, and the calibrated roll-flexure 
spring constant, Lọ, provide Lg, and I, from equations (30) and (31), 

Lo + aL from equation (32), and, for angle of attack other than zero, 
Lg from equation (33). 


(2) With the model free to yaw, the model is driven in roll at the 
yaw natural frequency. When the model has attained sufficient yawing 
amplitude, the drive motor providing the oscillatory rolling motion is 
stopped abruptly, resulting in a single-degree-of-freedom free-yaw oscil- 
lation. A time history of the decaying yawing motion is obtained by 
means Of a recording oscillograph. Again these data are obtained for 
both wind-off tunnel-evacuated and wind-on conditions, and with them and 
the calibrated yaw-flexure spring constant, Ny, equations (21) through 
(24) yield Iz, Mys, Ny - Ng, and Ng, respectively. 
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(3) With the model free to yaw, the apparatus is operated as a two- 
degree-of~freedom system at discrete frequencies throughout the operating 
range. At each frequency measurements of the roll position ©, yaw posi- 
tion y, roll-input torque T, and the phase angles € and A are obtained 
by means of the resolver circuits. Use of these data in conjunction with 
those obtained from (1) and (2) above provides ò + aNg and Ixyz from 
‘equations (25) and (26), Lg and Ly - Lg from equations (36) and (37). 


EVALUATION OF SYSTEM 


It should be noted that the aerodynamic moments produced by angular 
velocities are, in general, much smaller than the moments produced by 
angular displacements, and the disturbances inevitably present in any 
wind tunnel may be of the order of magnitude of the quantities to be 
measured. An analysis of the manner in which the effects of these dis- 
turbances are eliminated in the present system has been discussed in the 
section on instrumentation. The results of such analyses are not suf- 
ficient information to determine completely the over-all accuracy; however, 
a reliable indication of the accuracy can be obtained where it is possible 
to simulate the quantities to be measured with known inputs to the system. 
The results of tests of this type will be discussed. A further check of 
the system accuracy when operating under wind-tunnel conditions was obtained 
through the testing of a simple configuration for which, it was felt, cer- 
tain of the aerodynamic derivatives could be calculated fairly accurately. 
The results of these tests are discussed at the end of this section. 


Evaluation of the System With Simulated Aerodynamic Derivatives 


Cross derivatives No and Ly.- Tests with the simulated aerodynamic 
derivabliSl Were unde Ib Order te assess the accuracy of the system and 
the validity of the assumptions made in the linearization of the equations 
of motion. The aerodynamic cross derivatives > and Ly were simulated 
by means of a variable-speed gyroscopic motor as indicated in sketch (m). 
| The torque due to precession of a gyroscope is 

well known (see, e.g., ref. 8), and it is this 
property which was utilized in the present appli- 
cation. The gyroscope was mounted in place of 
y P the model, and, with the rotor oriented as shown, 
the simulated cross derivatives provided by the 
N gyro are given by 


| Ly = “Ng = 10 
Sketch (m) 
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where 


L rotor moment of inertia 
Q angular velocity of rotor 


The results of the gyroscope calibration are shown in figure 5 where 
the range of values given is expected to cover all possible wind-tunnel 
results. : It can be seen from figure 5 that the accuracy with which 
is determined is entirely satisfactory. It was found that the values 
of Ly obtained from equation (37) were unreliable, however, and these 
results are not presented. This was also found to be true for the calcu- 
lation of Lg from equation (36) and a re-examination of the method of 
determining these derivatives will be necessary. 


Damping-in-roll derivative ip. - The results of tests in which the 


aerodynamic roll damping I was simulated are shown in figure 6. This 
derivative was simulated in the manner indicated in sketch (n). A moment 
proportional to the angular velocity 
D was produced about the x axis 
by means of a magnetic damper acting 
through a flexure-pivot universal 
joint on an arm connected rigidly to 
the roll torque tube. The damping 


moment was generated by operating 
the system with a direct-current 
field excitation and a short-circuited 
armature. If the current generated Z 
in the armature and the armature 
resistance are designated as I and R, 
respectively, the simulated aerodynamic £ ` 
derivative may be expressed, for small 
angular displacement ©, as Sketch (n) 
‘ 2ncIFR 
Lọ = — ZZ 
Po W 


where ece is the conversion from joules to inch-pounds. With the damper 
inoperative, an equivalent damping factor was obtained for the internal 
friction in the roll flexure pivots, universal joint, and flexible mount- 
ing of the damper armature.. This quantity was subtracted from the damper- 
on data to give the results shown in figure 6. These results show the 
accuracy with which the system measures the simulated aerodynamic deriv- 
ative Lọ in the absence of tunnel disturbances. 


Damping in yaw and directional stability derivatives Cn, - Cno and 


Cng. No satisfactory method was found for introducing simulated values 


of aerodynamic yaw damping and directional stability into the system 
during the previously described tests, and it is, therefore, necessary 
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to analyze the results of the wind-tunnel tests, described in the follow- 
ing sections, to determine the accuracy with which these derivatives were 
obtained. The procedures for obtaining and reducing free-oscillation data 
are well known, as are the factors affecting the accuracy of these data. 
The uncertainties in each of the measured quantities present in the equa- 
tion defining Cnr - Cn3 and Cng have been calculated. The square root 
of the sum of the squares of these uncertainties indicated over-all 
uncertainties of the order of +0.055 and +0.057 for single values of 

Cn, - Cng and Cng, respectively. Not accounted for in these values, 
however, are possible effects of random air-stream disturbances and support 
vibrations which would appear as scatter in the data. The standard devia- 
tion of a number of observations at a given Mach number was, therefore, 
calculated; this was found to be of the order of +0.127 for Caa = Cg + 
and £0.028 for Cng- 


Evaluation of the System in the Wind Tunnel 


The tests described previously, in which the system accuracy was 
assessed from the standpoint of ability to measure simulated aerodynamic 
moments, were necessarily conducted under somewhat idealized conditions. 
It is, of course, impossible to simulate all of the loading conditions 
and disturbances to which the model is subjected in the wind tunnel, and 
these tests provide a measure of the accuracy of the system when there 
are essentially no extraneous disturbances present. It has been shown ` 
in the analysis of the previous section thet these disturbances, theoret- 
ically, cause no error; however, it is desirable to check the capabilities 
of the system operating in the wind-tunnel environment for which it was 
designed. 


An evaluation of the system in the wind tunnel necessarily depends 
on testing a model for which the aerodynamic derivatives may be calculated 
accurately. It is, of course, realized that an exact check on the system 
accuracy is not possible by this means since present theoretical knowledge 
is insufficient to permit exact calculations of the lateral dynamic sta- 
bility derivatives of fin-body combinations; a serious decrease of the 
system accuracy due to tunnel disturbances should, however, be evident. 


For the tests described below, a simple model was constructed for 
which, it was felt, the rolling-velocity derivatives Clo and Cnp could 
be estimated with a fair degree of accuracy. The model, which consisted 
of a vertical triangular fin mounted on a cylindrical body with an ogival 
nose (see fig. 7), was tested with the dynamic system in the Ames 6- by 
6-foot supersonic wind tunnel over a range of subsonic and supersonic 
Mach numbers at zero angle of attack. 


Estimation of stability derivatives.- Estimation of the stability 
derivatives was based on existing theoretical methods adapted to the 
particular model used in the wind-tunnel tests. The following is a brief 
description of the methods by which the derivatives were estimated. 
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In determining the roll-velocity derivatives, C7. and Cy, it was 
assumed that the fin acts as a portion of a rolling triangular wing as 
indicated in sketch (o). The solution for the resulting pressure dis- 
tribution is well known (see, e.g., ref. 11) and was used directly to 
calculate the damping-in-roll derivative C7 and the contribution of 
the fin to the cross derivative Cnp- Pressures induced on the body by the 
loading on the fin of course provide 
no rolling moment, but contribute ae 
directly to the yawing moment, and, Na Te 
therefore, must be considered in SO S 
estimating Cn. for the fin-body “w. T `s 
combination. Methods are available ISA s. 
for computing the interference load SSN 
distribution for fin-body combinations Se 
at supersonic speeds (refs. 12, 13, CAL 
and 14) but these are, in general, ‘\ l 
very laborious and for this reason N a 
the method described below was applied 
in the present case. Sketch (o) 


It was assumed that the pressures induced on the body could be cal- 
culated as though the body were a flat plate at zero angle of attack 


adjacent to the fin. Reference 15 _ bz 
gives the solution for a case similar Q a vV 


to this, namely the pressure distribu- 
tion induced on a flat plate at zero 
angle of attack adjacent to a tri- 
angular sector at a constant angle 

of attack. The angle-of-attack dis- 
tribution on the fin may be repre- 
sented as the sum of a constant angle 


of attack (= =) and one which 


V 
increäses linearly in the spanwise 


direction (a = Be) as shown in 
sketch (p). — Z 


The pressure induced on the 
body by the constant angle-of-attack p(z-r) 
distribution is given directly by the a Qz E aa 
solution in reference 15. The induced 
pressure resulting from the effective 
twist due to rolling velocity may be 
found by integrating across the span 
of the wing. This integration cor- 
responds to the superposition of an 
infinite number of the lifting sectors : Sketch (p) 
along the span, each sector having an infinitesimal angle of attack. The 
yawing moment due to the induced pressures on the body was added to that 
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due to the loading on the fin to obtain the total yawing moment due to 
roll velocity. Because of the approximations necessary in the above 
estimate, the yawing moment due to rolling velocity was also calculated 
using the slender-body theory of reference 16 in order to provide a 
check of the order of magnitude. 


In estimating the yawing derivatives it was assumed (a) that, in 
lieu of a more exact theory, the yawing moment of the body is that given 
by slender-body theory which, of course, is known to be incorrect at the 
higher supersonic Mach numbers and (b) that the yawing moment of the fin 
due to yawing velocity is that given by the moment, about the appropriate 
axis, of the loading on a portion of a pitching wing. Sketch (q) shows 
the portion of the wing referred to as a cross-hatched area and the 

pitching axis as the axis of yaw. The load- 
ct ing was calculated according to reference 1f 
for the arrangement as shown to determine 


NS me 1 the yawing moment due to yawing velocity. A 
Sg 5S similar method was used for estimating the 
S PA yawing moment due to sideslip angle. The 


ee We Slender-body-theory results were obtained 
“ys from reference 18. The above methods, of 


` course, applied only to the derivatives for 


t NS y supersonic Mach numbers. No methods were 
| NOT available of sufficient accuracy to permit 
\! the estimation of the derivatives at subsonic 
v i 
Sketch (q) Mach numbers. 


Comparison of estimated stability derivatives with measured valués.- 
During the wind-tunnel tests the rolling derivatives Clp and Cnp were 


obtained over a range of frequencies from Į to 20 cycles per second. No 
significant variation of these quantities with frequency was found, and 
the values which are presented in figures 8 and 9 are the arithmetic 

mean of the values obtained. The values of Cn, ~ Cng and Cna are those 
which were obtained at a single frequency through free-oscillation tests. 
A comparison of the estimated and measured rolling derivatives is given 
in figures 8 and 9 for the supersonic Mach numbers investigated. The 
measured values are also shown for the subsonic Mach numbers. The agree- 
ment between the experimental and estimated values indicates that there 
were no serious effects of wind-tunnel disturbances on the system accuracy. 
Figure 10 shows the measured and estimated damping in yaw. Here the 
agreement is poor, particularly at the highest Mach number. The same is 
true of the directional stability parameter Cne, shown in figure ll, 
Where the agreement is satisfactory except at the highest Mach number. 
The methods used in estimating the theoretical values of the derivatives 
may not be reliable and have not been adequately checked experimentally. 
A large part of the discrepancy between theory and experiment is believed 
to be due to deficiencies in the theory. 
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CONCLUSIONS 


A study of a system designed to measure the static and dynamic 
lateral stability derivatives has lead to the following conclusions: 


l. Tests in which aerodynamic derivatives were simulated have shown 
that the apparatus correctly measures known inputs of the order of magni- 
tude anticipated for airplane models. 


2. The investigation has shown that the aerodynamic derivatives 
associated with the primary mode of operation of the system, that is the 
rolling derivatives, and the directionel-stebility and damping-in-yaw 
parameters may be obtained with satisfactory accuracy. The data-reduction 
equations from which the rolling moment due to yawing velocity and due 
to sideslip angle were determined yielded results which were unreliable 
and prevented evaluation of Cis - Uig and Gig! 


3. Tests of the apparatus utilizing a simple model showed that the 
wind-tunnel environment apparently had no deleterious effects on the 
system accuracy. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Sept. 24, 195h. 
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Figure 1.- Illustration of coordinate system. 


Displacements and moments are shown in the positive 
sense. 
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(a) Schematic drawing. 
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Figure 2.- Illustration of dynamic testing apparatus. 
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Figure 2.~ Concluded. 
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Figure 3.- Effect of simulated lift load on yaw flexure pivot spring constant. 
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Figure 4.+ Effect of simulated lift load on the internal damping of the yaw flexure pivots. 
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Figure 5.- Comparison of gyroscope input and measured output of the simulated aerodynamic cross 
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Figure 6.- Comparison of magnetic-damper input and measured output of the simulated aerodynamic 
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Figure {.- Model used for wind-tunnel tests, 
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Figure 8.- Variation with Mach number of experimental and estimated values of damping-in-roll 
coefficient, 
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Figure 9.- Variation with Mach number of experimental and estimated values of cross-derivative 
coefficient, Cp» 
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Figure 10.- Variation with Mach number of experimental and estimated values of damping~in-yaw 
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Figure ll.- Variation with Mach number of experimental and estimated values of directional atability 


NACA-Langley - 1232-3-54 - 1000 


